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We present field theory simulations of a model with Z2 x [/(I) symmetry in (2+l)-dimensions. This 
model has two discrete vacua, allowing for domain walls, and also a conserved Noether charge. For 
initial conditions in which the field is placed randomly in one of the vacua and given a homogeneous 
background charge, we find that the number of walls does not scale in the standard way. We argue 
that the Noether charge and current become localized on the walls, forming kinky vortons, (2+1)- 
dimensional analogues of cosmic vortons. These loops of wall can be long lived, or even stable, 
depending on their precise characteristics. We suggest that our simulations illustrate a possible 
mechanism for dynamical frustration of domain wall networks and that cosmic vortons will form 
naturally in (7(1) x (7(1) models. 



I. INTRODUCTION 

Domain walls are field theoretic solutions which interpolate between discrete vacua. They can be of interest in 
the context of particle physics and cosmology where the field is relativistic and fundamental [1], or in condensed 
matter where the field is non-relativistic and represents an order parameter of the system [2] . Once formed at a phase 
transition, the number of walls is expected to scale cx in relativistic systems [3 HI EJ El [7] which will be the focus of 
the present paper; essentially a network will lose energy as quickly as possible within the bounds imposed by causality 
and the causal horizon is the only length-scale in the problem. Similar standard results exist for non-relativistic 
systems |H]. 

Cosmic vortons are loops of superconducting string which are stabilized by the existence of charge and current 
which form a condensate coupled to the vortex field [TUj- Recent numerical work has investigated the stability of such 
configurations in (3+l)-dimensions |TTJ [T31 E] , but due to the large variety of length-scales in the problem, this 
has only allowed quantitative validation of the standard picture for their formation and evolution in a small range of 
parameter space. 

Kinky vortons are (2-t-l)-dimensional analogues of cosmic vortons which replace the superconducting vortex with 
a superconducting domain wall {15j (or kink) . In this lower dimensional case extensive numerical investigations have 
been possible. In addition due to the fact that exact analytic solutions are possible for the superconducting condensate, 
independent of the charge and current, this has allowed analytic investigation of their stability properties in the thin 
ring limit |[16j developed by Carter and collaborators [171 UHl HH EHl Ell l22l ES] . This work appears to suggest that 
the basic tenets of cosmic vorton evolution, that they might eventually come to dominate the energy density of the 
Universe, are correct (see, for example, refs. [24l [25l [26] ). 
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In this paper we will investigate the evolution of networks of these superconducting domain walls formed from 
random initial conditions coupled to a homogeneous background charge. In particular we will investigate the scaling 
properties of the number of walls as a function of time. We find that the effects of the conserved charge and current, 
which localize on the walls, appear to slow down the natural propensity of kinky-loops to collapse under their own 
tension, violating the standard law of scaling. Some of the loops appear to have similar properties to stable kinky 
vortons suggesting that such objects can form rather naturally and, although we cannot confirm their absolute stability 
due to our limited dynamical range, they are, at the very least, long-lived. 

The scaling properties of domain wall networks have another interesting cosmological application. It has been 
suggested that they may form a frustrated lattice which could be a candidate for a dark energy component with 
w — P/ p — —2/3 |27l |28]. However, all simulations of domain wall networks, even those with junctions (for example 
refs. Ill |5j [6l [7j), appear to be compatible with the standard scaling law - although they all have relatively poor 
dynamical range compared to cosmological scales. In a separate paper |29j we will investigate how similar models to 
those considered here, which also allow for junctions, behave in an analogous way, but with some extra complications. 



II. MODEL AND NUMERICAL IMPLEMENTATION 



The (2-t-l)-dimensional model we will investigate is described by the Lagrangian density 

C = d^c^d,^ + d^ad.-a - ^ (0^ - r^)' ~ ^ {\-? - €)" " /30Vl^ (1) 

which has a global Z2 x t/(l) symmetry. Throughout the subsequent discussion we fix the model parameters — 
Act = 2,/3 = 1; see ref. [15 for a discussion of parameter values in this theory, and the motivation for this particular 
choice. The real scalar field x) is the broken domain wall- forming field and the complex scalar field cr(t, x) will 
act as the superconducting condensate. 

The equations of motion admit a y-directed domain wall solution |15j 

= ?70tanh {rn^x\J\ - a - x) ■, (2a) 
a = |(7|e*('=^+"*\ (2b) 



\a\ = \/2"?70Wa- ^ + x sech^^m^xyr^^^i^^) , (2c) 



where we have defined a = {rj^/rj^) , x = x/"''0' 'X^ — ^'^ ~ 

2 ''^<t> 2 2 '^c 2 I 

Existence of a solution requires that 1/2 — a<x< 1 — a and this solution explicitly exhibits many of the expected 
properties of superconducting solitons such as current quenching. 

One can use this solution to construct ring solutions if the radius, R, is sufficiently large that corrections due to 
the curvature are negligible. Such kinky vorton solutions can be specified by the winding number N with k = N/R, 
and the conserved charge of the condensate field is given by 

Q ^ J d'^xpQ 

The current is given by 



J = 1 (ctVct - ctVct) . (4) 



The stable radius Ri, can be computed for sufficiently large values of N and Q 



3 



We evolve the equations of motion in (2+l)-dimensions by discretizing on a regular square grid of points and 
grid-spacing Ax. The time evolution is performed using a second-order leapfrog algorithm, with time-step size At and 
spatial derivatives are discretized to fourth order. Typically we will use P = 4096 and Aa; = 0.5, At = 0.1, although 
some simulations have used P — 1024 in order to reduce computational requirements. All simulations presented here 
use = i,rjcr = \/3"/2 (such that a — 3/4), but similar results are found for all values of a compatible with the 
existence of a kinky vorton solution. We use periodic boundary conditions, which sets a light-crossing time: the 
simulation time taken for an emitted signal to interact with itself after passing through the periodic boundaries, given 
by T = iPAx. Our results on the scaling of the number of walls are only valid upto this light-crossing time; however, 
structures that persist after this time are still of interest, particularly when they appear to be stable. 

Initial conditions are chosen in order to represent a phase transition in the (/)-field and a condensate with a uniform 
charge density in the cr-field. We note that it is not our intention to model a realistic phase transition by which kinky 
vortons might be formed, rather to do something which is random and relatively painless from a practical point of 
view. At each grid-point the 0-field is placed randomly into one of the two points on the vacuum manifold, with 
= 0. The cr-field is set to ct = Ae"^* where A and uj are real constants for which the initial charge density is 
Pq(0) — A^oj. All simulations presented here have lo — 1.0, and A is varied to change the initial charge density; we 
have checked that similar results are found if we fix A and vary lj. Although these initial conditions are random, 
they will contain unphysical energy gradients between adjacent grid-points. In order to ameliorate this problem we 
introduce a damping term which only acts on the (/)-field for the first 200 time-steps in our simulation. This results 
in a reduced energy and aids condensation into smooth domains, without reducing the total charge which remains 
approximately constant. 

In order to study the scaling properties we need a method which can be used to estimate the number of walls in 
the simulation from the fields. We do this by first deciding which of the two minima (p is nearest at each grid point. 
If either of the adjacent grid points (only above and to the right to avoid duplication) return different minima then 
the wall count is incremented. 

III. RESULTS 

We have peformed simulations with P = 4096 upto t — t for four different values of A which correspond to initial 
values of /Oq(0) = 0.25, 0.09, 0.01 and 0. The evolution of the fields are presented in Fig.jljat times which are multiples 
of 2 from t = 80 upto 1280 (which is just beyond t = t). If the network is scaling, then one would expect the size of 
the domains to increase by a factor of 2 from one plot to the next. The first thing to note is that the zero and low 
charge case (pq(0) = 0.01) look very similar, whereas the higher charge cases (pq(0) = 0.09 and Pq(0) = 0.25) are 
qualitatively different, with much smaller domains particularly in the case of Pq(0) = 0.25. The number of walls as 
a function of time are presented in Fig. [2[ This confirms our observations about Fig. [lj for pq(0) = and 0.01 the 
number of walls appears to be oc t~^, whereas for pq(0) = 0.09 and 0.25 it is clearly oc t~'^ where the scaling exponent 
7 <C 1. In fact for /Oq(0) = 0.25 the number of walls appears to freeze to a constant value at late times. 

This is the main result of this paper. We find that domain walls formed in the presence of charge do not scale in 
the standard way and it appears that, if the initial charge density is sufficiently large, the configuration appears to 
freeze in. This is further illustrated by Fig. |3]which presents simulations with /Cq(0) = 0.25 and pq(0) = 0.09 beyond 
t = T for P = 1024, which was used in order for the simulations to be completed in an acceptable period of time. 
We see that, although the simulations contain a lower number of domains due to the smaller grid, there is very little 
evolution of the structures beyond t — 2t. 

In order to investigate the mechanism by which this non-standard scaling behaviour is achieved we have presented 
a number of different quantities for the simulation with /5q(0) — 0.09 at t = 640 in Fig. |4] We note that in each 
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FIG. 1: Evolution of the (^-field for Pq(0) = 0.25, 0.09, 0.01, (from top to bottom). The colours correspond to each of the two 
points on the vacuum manifold (red/light: 4> = +1; blue/dark: cj> = —1). The images are at f = 80,160,320,640,1280 (left 
to right) and we have used P = 4096. There is a clearly visible difference in the long-term behaviour of the low charge cases 
Pq(0) = and 0.01, and the higher charge cases pq(0) = 0.09 and 0.25. 

case we have tuned the colour scheme to highlight the important features of the plots. First in (a) we present the 
potential energy density, ppE. This is seen to be localized around regions corresponding to where (j) goes from positive 
to negative confirming that domain walls have been formed. There is one very large loop which goes through a 
number of the periodic boundaries. Since it is possible to draw lines top-to-bottom and left-to-right which cross the 
boundaries but do not cross any of the walls, this loop does not wrap around the toroidal geometry imposed by the 
periodic boundaries. There are also a number of other smaller loops, most of which are clearly close to circular, and 
one which is in the middle of the plot which has a particularly non-regular shape. 

In (b), (c) and (d) we have presented 3?((t), pq and the modulus of the current |Jp. All three quantities are mainly 
localized around the wall structures seen in (a). We see that ^{a) alternates from positive to negative as one traverses 
a particular loop, indicating that the cr-field is twisting along the wall; reminiscent of that seen in a kinky vorton. 
The fact that the charge and current are localized on the wall reinforces this conclusion. 
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FIG. 2: Evolution of the number of domain walls for various Pq(0) with P — 4096. For Pq(0) = and 0.01 we see that the 
number of walls appears to scale oc , whereas the evolution is modified for large values of pq(0). 




FIG. 3: Evolution of the <^-field for pq(0) — 0.25 (top), 0.09 (bottom), at t = r, 2r, 3r, 4r (from left to right). Due to dynamical 
range constraints we have used P — 1024 for these simulations. We see that there is very little evolution beyond the light 
crossing time, suggesting that the structures have frozen in. 

We are, therefore, led to the conclusion that it is the charge and current on the walls which are counteracting the 
string tension in the same way as in kinky and cosmic vortons, and preventing the network from losing energy as fast 
as causality will allow as is necessary for the standard scaling regime. In fact, as we will discuss below, we believe 
that only two of the loops (the very large loop and the smaller unusually shaped loop in the middle of the plot) are 
close to stable, but the others are much longer-lived than they would be in the absence of charge. 

In order to investigate whether the long-lived loops are indeed kinky vortons we have attempted to estimate the 
values of N, Q and R for some of the loops in the pq{0) = 0.09 simulation. The loops we have selected are illustrated 
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FIG. 4: Results of the simulation with Pq(0) — 0.09 at t = 640 : (a) potential energy density, ppe using colouring routine where 
purple is the maximum value, grey is where less than 10% of its maximal value; (b) K(cr) where red/light and blue/dark denotes 
K(a) positive/negative, and grey where K(it) < 40% maximum value; (c) charge density, pq, where red denotes positive, grey 
denotes pq < 10% maximal value; (d) the modulus of the current | with blue denoting | J|^, grey denotes jJp < 2% maximal 
value. All images are coloured on a gradient until the specified threshold. Counting loops reveals one large loop and a number 
of other smaller loops none of which wrap around the toroidal geometry. It is clear that ppB, pq and | Jp are localized on the 
domain walls, and that K(a') oscillates between positive and negative values as one traverses a loop of wall. 
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FIG. 5: Properties of three of the loops in the Pq(0) — 0.09 simulation. TV is found by counting the number of red-blue pairs 
along a loop; Q is computed by integrating the charge density over a region 10 grid-squares from the loop; and R is computed 
from the area of the loop, i?* is the radius of a stable kinky vorton with the given and Q. The highlighted loops are those 
for which the table gives the properties. The loops are labelled from top left to bottom right (1-3). The fraction of the charge 
which resides on all loops to the total charge in the box (i.e. not just those highlighted) is 7% and 6% for the two time-steps 
shown. 



in Fig. |5] They have been chosen since they are localized and present in the simulation when t = 640 and t = 1280. 
Two are close to being circular and the other is the unusually shaped loop already discussed above. The radius is 
estimated by calculating the area enclosed by the loop and equating this with ttR'^; obviously this is only a crude 
approximation. Using the area is more accurate than the circumference since this is susceptible to errors induced by 
small-scale structure. The value of Q is computed by including charge which is located within 10 grid squares from a 
point where a wall is identified. The value of N is estimated by eye from Fig. Qb), by counting the number of red to 
blue transitions. The values are presented in the table in Fig. [Sj along with the stable radius, i?^,, computed for the 
estimated values of Q and N, using the formula given in ref. [T5] . 

We see that for loop 2 the values of R and i?* are consistent with each other (within the somewhat rough calculation 
we have used for such a non-circular loop), and that the size of the loop reduces only a small amount with time. This 
indicates to us that this is close to a kinky vorton. For the other two loops (1 and 3) the estimated values of R are 
somewhat larger than i?^, and the values of N have reduced quite a bit between the two timesteps. We conclude that 
the loops are not likely to stabilize, but that their collapse is being inhibited by the charge and current on the loops; 
however, one should note that the radii they are collapsing from is greater than the stable radii - perhaps indicating 
that the loops could attain the stable radii, given sufficient simulation-time. We note that all the loops formed are in 
the electric regime, that is Q ^ 47r7V. 

It is difficult to reproduce a similar analysis for the case of Pq(0) = 0.25. Although there are actually many more 
loops, their density is much higher and they interact with each other, making it difficult to identify loops which 
are present at both timesteps. Moreover the equivalents of Figs. |4|^b),(c) and (d) are much more confusing and, in 
particular, it is difficult to estimate N with any confidence for any loops. We have already pointed out that in this 
higher charge simulation the situation is much closer to that of dynamical frustration with the high density of walls 
leading to a quasi-stable network. However, true frustration will probably require the underlying field theoretic model 
to include junctions. 
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IV. CONCLUSIONS 

Our work has shown that the evolution of superconducting domain walls can lead to a network which does not 
necessarily scale in the standard way. In the case of a uniform background of charge as discussed in this paper, there 
appears to be a critical value of pq{0) above which the scaling exponent 7 < 1 and possibly a higher charge where 
this becomes zero. 

We have conclusively shown that charge and current are localized onto the loops and that they can be long-lived. 
Some appear to be close to stable and have the properties of kinky vortons. In the case of P = 4096 and Pq{0) = 0.09 
it was possible to synthesize just two loops which were close to stability and a host of others whose collapse was 
inhibited by the existence of current and charge. For lower values of P, for example 1024, there is only one large loop 
which is left at the end. We would presume that more stable loops could be found for larger values of P and these 
small numbers are just a reflection of the low dynamic range. 

Although kinky vorton solutions have been shown to exist and cosmic vortons have also been invesitgated in a 
narrow parameter range, their creation in a truely random scenario might have seemed a little implausible since they 
are dynamically, rather than topologically, stable. However, our conclusions appear to suggest that the formation 
of vortons as part of a more realistic scenario is at least possible. Of course, as we have seen most loops will be 
unstable and will ultimately collapse, but the cosmic vorton scenarios only require a very small number of loops to be 
absolutely stable. Unfortunately, to investigate the distribution of vorton sizes formed in a realistic scenario requires 
much greater dynamical range than that used in the current study. 

Our simulations also appear to show possible direction for those interested in forming a frozen lattice from domain 
wall networks. For example, they form a caveat to the claims made in ref. [B] which state that dynamical frustration 
is impossible. Nonetheless, much work is necessary to turn this embryonic idea into a full theory of dark energy. 
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